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Abstract The observational basis for asteroseismology 
is being dramatically strengthened, through more than 
two years of data from the CoRoT satellite, the flood 
of data coming from the Kepler mission and, in the 
slightly longer term, from dedicated ground-based fa- 
cilities. Our ability to utilize these data depends on 
further development of techniques for basic data anal- 
ysis, as well as on an improved understanding of the 
relation between the observed frequencies and the un- 
derlying properties of the stars. Also, stellar modelling 
must be further developed, to match the increasing di- 
agnostic potential of the data. Here we discuss some 
aspects of data interpretation and modelling, focussing 
on the important case of stars with solar-like oscilla- 
tions. 

Keywords asteroseismology; stellar internal physics; 
stellar structure; stellar evolution; helioseismology; so- 
lar composition 



1 Introduction 

In the last two decades data have become available 
on solar oscillations, which have allowed detailed infer- 
ences to be made on the prop erties of the solar interior 
( Christensen-Dalsgaard 20021 ). This has led to the de- 
termination of the solar sound speed with considerable 
accuracy, detailed tests of the equation of state of the 
solar interior and a determination of the solar internal 
rotation, the origin of which is still poorly understood. 
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The inferences about the solar internal structure obvi- 
ously provide a strong test of modelling of stellar evolu- 
tion. Although earlier solar models were in reasonable, 
but far from pe rfect, agreement w ith the helioscismic 



determinations ( Gough et al. 19961 ). recent redetermi 



nations of the solar surface abundance have led to a 
very substantial increase in the discrepan cies between 
the Sun and the model ( for reviews, see lOuzikl 120061 . 
20081 iBasu fc AntialEool see also Section E]). 



The Sun is only one star at a specific stage of its evo- 
lution, and its structure is relatively simple, compared 
with other stars. Thus it is of obvious importance to 
extend seismic investigations to other stars. This is be- 
coming possible, thanks to the development of ground- 
based instrumentation and the launch of three space 
missions that are providing extensive data on stellar 
oscillations. Thus it is important to consider the tools 
available for the analysis of such asteroseismic data and 
summarize the results that have already been obtained. 

Here we concentrate on solar-like oscillations; these 
are characterized by being intrinsically damped and 
excited stochastically by near-su rface convection (e.g., 
Houdek et alll999tlHoudekll2006h . Thus they arc found 
in relatively cool stars, with significant outer convec- 
tion zones, both on the main sequence and on the 
red-giant branch, as rec ently confirmed specta cularly 



by the CoRoT mission (|De Ridder et al.ll2009l ). The 



modes relevant to observations of distant stars have 
low spherical-harmonic degree and hence are generally 
high-order acoustic modes; however, in evolved stars 
they may take on a mixed character, with a component 
corresponding to an internal gravity mode in the deep 
interior of the star, v ery substantially enhancing t heir 
diagnostic value (e.g. 
also Section [ 



Christensen-D alsgaard 2004 see 



At the simplest level, the oscillation frequencies are 
determined by the global properties of the star, more 
specifically the mean density p oc M/R 3 , where M is 



2 



the mass and R the surface radius of the star. When 
combined with additional, nonscismic, data and stel- 
lar modelling, the stellar radius can be determined 
with substantial p recision and, one may hope, accuracy 



( Stello et al.ll2009l ). Such a determination of the radius 
is particularly important for observations of extra-solar 
planets using the transit technique, where determina- 
tion of the radiu s of the planet requi res knowledge of 
the stellar radius (jKjeldsen et al 1 l2009h . From the point 



of view of stellar modelling, however, it is of substan- 
tially greater interest to obtain information about the 
stellar interior. From low-degree acoustic modes a mea- 
sure can be obtained that is sensitive to the properties 
of the stellar core and hence, with additional assump- 
tions, provides a measure of the amount of hydrogen 
consumed during stellar evolution and thus of the stel- 
lar age. With sufficiently detailed and accurate obser- 
vations additional information can be obtained about 
the stellar interior, including inferences of the sound- 



star (e.g., 


Gough & Kosovichev 


1993; 


Basu et al. 


2002 


Roxburgh 2004). 



An obvious requirement for asteroseismic inferences 
is the availability of good observational data. This is 
increasingly being met (see Section EJ) . Ground-based 
observations of Doppler velocity have been possible for 
a few stars, at level of sensitivity of a few cms -1 . 
Space-based asteroseismic observations were started 



serendipitously by the WIRE satellite (jBuzasi et al 



20051 lBrunttll2007[ ). followed by the launch in 2003 of 
the highly successful Canadian micr o-satellite MOST 

The French-led 
20091 ). launched 



( Walker et all 2003; Matthews 2007). 



CoRoT satellite mission ( Baglin et al 



at the end of 2006, has provided photometric oscillation 
data on a large number of stars, and high-quality data 
on an even more extensive stellar sample is promised 
by the NASA Keple r mission launched in March 2009 



(jBorucki et al.l l2009h . Further observational projects 



promise a substantial increase in the quality of data 
available for asteroseismic analysis. However, the inter- 
pretation also depends on an adequate understanding 
of the relation between the observed oscillation prop- 
erties, in particular the frequencies, and the underly- 
ing properties of the star, in order to design optimal 
diagnostics based on the observations. Furthermore, 
the physical interpretation of the results relies on stel- 
lar modelling, relating the physics of stellar interiors to 
stellar structure. Thus a crucial aspect of the analysis is 
the reliable computation of stellar models. This should 
obviously ensure numerical accuracy; an important step 
towards this goal is the detailed comparison of stellar 
models and oscillation frequencies computed by inde- 
pendent codes that has been undertaken in the ESTA 



collaboration, initially a s part of the CoR oT project 
(jLebreton et all l2008al[b1 ; iMova et al1l2008h . However, 
an equally important aspect is to ensure that the un- 
avoidable approximations in the modelling are imple- 
mented in a consistent and accurate manner; only in 
this case can the asteroseismic analysis provide a reli- 
able assessment of these approximations and, ideally, 
point to ways to improve them. Finally, the results 
should be presented in a statistically well-defined man- 
ner, such that the significance of the inferred properties, 
and of any deviations from the assumed models, can be 
assessed. It is probably fair to say that we are still some 
way from meeting these requirements. 

In the present paper we discuss some aspects of the 
diagnostics of solar-like oscillations, based on asymp- 
totic relations and computed properties; this in par- 
ticular emphasizes the important probing of stellar 
cores. In addition, we consider some examples of un- 
certainties in stellar internal physics which have been, 
or may be, amenable to seismic inferences. More 
extensive overviews of helio- and aster oseismi c diag - 
nost ics were provided by, f o r exa mple, iGoughl (|1993l ) 
and IChristensen-Dalsgaardl (|2004f ). and a comprehen- 
sive discuss i on of asteroseismology can be found in 
Aerts et al.l <|2009h . 



2 Asteroseismic diagnostics 

2.1 Simple asymptotics of acoustic modes 

The observed oscillation modes in main-sequence solar- 
like stars are driven predominantly by the vigorous tur- 
bulence in the superficial stellar layers. This results in 
a rich acoustic pulsation spectrum with a characteris- 
tic comb-like frequency structure. In distant stars only 
low-degree modes can be observed; however, the ob- 
served modes are generally of high order allowing us to 
use asymptotic theory. For modes satisfying l/n — » 0, 
where n is radial orde r and I degree, the frequencies v n i 
can be estimated as (|Vandakurovlll967l ; lTassoullll980l ; 
Goughlll986h : 



/ 



v nl ~ I n + - + e i/o - 



AL 2 -B 



where 



i/q 



21 * 

o c 



^o+OK 4 ), (1) 



(2) 



is the inverse of twice the sound travel time between 
the centre and acoustic surface, c is the sound speed, 
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and 



A 



4ir 2 v 



c(R) 
R 



dc dr 
dr r 



(3) 



The frequency-dependent coefficient e is determined by 
the reflection properties of the s urface layers, as is 
the small correction term B (e.g., iGoughl 119861 ). and 
L = sfl{l + i). The value of vq can be estimated from 
taking the average (over n and I) of the so-called large 
frequency separation 



l = V n+ \l — V n i 



(4) 



between modes of like degree and consecutive order. 
The resulting comb-like frequency structure consists of 
modes of odd degree falling approximately halfway be- 
tween modes of even degree. The last two terms on 
the right-hand side of equation (JTJ) lift the degeneracy 
between modes with the same value of n + 1/2, lead- 
ing to the so-called small frequency separation 5v n i = 
v n i — The small frequency separation mea- 

sures the extent to which the integrated history of nu- 
clear transmutations has modified the structure of the 
energy-generating core. The mean small frequency sep- 
aration, averaged over n (indicated by angular brackets) 
becomes in the limit l/n — > 



(SVnl) = (v n i - Vn-ll+2) - 2A(2l + 3) 



Wnl) ' 



(5) 



and is therefore proportional to the coefficient A 
(Eq. [5]) . The coefficient A measures predominantly the 
stratification of the energy-generating core through the 
gradient of the sound speed and hence is sensitive to the 
chemical composition ther e and consequently is an indi- 
cator for stellar age (e.g ., Christensen-Dalsgaardl f 984 , 
199, 4 lGoughlll986l . l200ll : IUlrichlll986l : lHoudek fc Cough 
2008h . 

From Eq. ^ it follows that, to this asymptotic order, 



D„ 



Vnl — Vn-\ 1+2 

21 + 3 



(0) 



is a function of frequency. The same is true for the 
following ratio: 



£,(1) _ \Vn-l I + V n i)j2 — V n -H+i 



(7) 



which is also a useful diagnostic, particularly when only 
modes of degree and 1 can be identified from obser- 
vations. 
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Fig. 1 The internal phase shift 81 (cf. Eq. [8]) for a model of 
the present Sun, as a function of cyclic frequency. Continu- 
ous curve: I = 0; dotted curve: I — 1; dashed curve: I — 2; 
and dot-dashed curve: I — 3. 

2.2 The inner phase shift 

Further insight into the diagnostic potential of low- 
degree solar-like oscillations has been developed by 
Roxburgh and Vorontsov. They represent the effect 
of the structure of the core by an internal phaseshift 
Si(uj) (jRoxburgh fc Vorontsovlll994l ), defined such that 
the asymptotic behaviour of the acoustic modes can be 
expressed as 



0*0 



1/2 



-p 



Ap sin (ujt 



where p is density, p' is the Eulcrian perturbation to 
pressure, lu = 2-kv is angular frequency and 



r dr 



(9) 



is the acoustic distance from the centre, corresponding 
to the radius r. They noted that 5i can be computed as 
a function of frequency by matching partial solutions to 
the oscillation equations to the asymptotic expression in 
Eq. (JSj) . We have implemented this by considering solu- 
tions satisfying both the central and surface boundary 
conditions, but allowing a discontinuity in the solution 
at a suitably chosen point near the surface. Figure [1] 
shows the results for a model of the present Sun. 

From Eq. §E§ and a similar asymptotic representation 
of the solution in the outer parts of the model, the 
oscillation frequencies can be obtained as 



Vrd = V a 



n+-+7T 1 {5t(uj nl ) - 8l{fjJ n l)) 



(10) 



where a(u>) is a surface phase function (se e also 
Christensen-Dalsgaard fc Perez HernandeJl992h which, 
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Fig. 2 Dimensionless squared sound speed, in units of 
GM/R, as a function of fractional radius, for models of 
a Cen A (solid) and B (dashed) 



for low-degree modes, is independent of degree. This 
expression is clearly closely related to the simple asymp- 
totic expression in Eq. |T]) above, with the small fre- 
quency separation being related to 5i through 



v n l — V n -H + 2 — —{Si+2 — Si) 

7T 



(11) 



Roxburgh fc Vorontsovl (|1994l ) noted that, in the 



simplest asymptotic approximation, Si can be expressed 
as 



Si(uj) = 8 (uj) + 1(1 + 1)S{uj) 
where 

2uj 



'c(R) 


f R dc dr" 


R 


Jo dr r 



(12) 



(13) 



Given Eq. (fTTjl . this evidently leads to Eq. Q for the 
average small separation, with A given by Eq. (|3"j). 

The simple asymptotic behaviour in Eqs |T]) and (|3|) , 
or Eq. (fT2|) and Eq. (fT3|) . assumes that the equilibrium 
structure is varying slowly with position. As a star 
evolves on the main sequence, the increase in the cen- 
tral mean molecular weight leads to a local decrease in 
the sound speed, as illustrated in Fig. [5]for the compo- 
nents of a Cen binary system. At the age of 6.98 Gyr of 
the system the B component, with a mass of 0.928 M Q , 
is relatively unevolved, while the A component, with 
a mass of 1.111 M©, is nearing the end of central hy- 
drogen burning. The effect is evident in the difference 
in the behaviour of the dimensionless sound speed and 
is reflected in the scaled small separations shown in 
Fig. [31 While the scaled separations, in accordance with 
Eq. H}, are essentially just a function of frequency for 



_ 6 

4 . 



#: I = £:l = l ■: I = 2 
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Fig. 3 Scaled small frequency separations (cf. Eqs [6] and 
0), for models of a Cen A (panel a) and B (panel b). Filled 
symbols: , for I = 0, 1 and 2; the degree is indicated by 
the symbol, as shown in the figure. Open symbols: D n i, for 
/ = and 1. Note that according to Eq. |I} (see also Eq.O 
these quantities are expected to be functions of frequency 
alone 



a Cen B, for the A component there is very substan- 
tial scatter, reflecting the rapid variation of the sound 
speed in the core. 

This behaviour is illustrated in Fig. [4J in terms of 
the scaled differences [1(1 + 1)] -1 (<5/ — So) in the internal 
phase shift, which according to Eq. (fT2"]l is expected 
to be a function of frequency. This is clearly approxi- 
mately satisfied for the B component, whereas for the 
A component the results for I = 1 deviate strongly 
from the expected asymptotic behaviour. Interestingly, 
the differences for I = 1 show a distinct oscillatory be- 
haviour whic h is related to th e base of the convection 
zone (see also Roxburgh 20091 and Section 12.31). 

It was noted by Roxburgh fc Vorontsov (|2003 ) that 
the uncertainty in the modelling of the near-surface lay- 
ers (see Section |2~51) has a significant effect on the small 
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Fig. 4 Scaled differences in internal phase shift, for models 
of a Cen A (top) and B (bottom). Dotted curve: I = 1; 
dashed curve: I — 2; dot-dashed curve: I = 3. Note that 
according to Eq. (|12[1 these differences are expected to be 
functions of frequency alone 



frequency separation. They pointed out that this effect 
is strongly suppressed by considering instead ratios such 



VnQ — Vn-\ 2 
V n -1 1 — v nl 
Vn\ — V n -\ 3 



(14) 



They also demonstrated that to this accuracy the ratios 
are related to the internal phase shifts by 



r„o = {S 2 - 5 Q )/n , r nl = (S 3 - 5i)/n 



(15) 



This is illustrated in Fig. [S] for a model of a Cen 
A. The diagnostic potential of these separation ratios, 
and their insensitivity to even very substantial diffcr- 
ences in the surfa c e lay ers, was further analy zed by 
Oti' Floranes et all |2005h and iRoxburgbl ^OOol ). 
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Fig. 5 Frequency-separation ratios (cf. Eq. 1141 symbols) 
and difference (<5;+2 — 8i)/n in internal phase shifts (curve), 
for a model of a Cen A. Top: / = 0, bottom: / = 1. 



One may hope that further analysis of the internal 
phase function can provide insight on the diagnostic 
potential of the fre quency separations, beyo n d the sim- 



plest asymptotics. iRoxburgh fc Vorontsovl (jl994f ) de- 
veloped relatively simple expressions for Si, in terms 
of the structure of the core, although these have ap- 
parently not been exploited in de tail so far (see also 
Roxburgh fc Vorontsovl[2000ll2007h . 

As discussed in Section 14.21 the diagnostics of con- 
vective cores is particularly important. The mixing 
in the core leads to sharp variations in composition 
and hence sound speed there, causing strong depar- 
tures from the simple asymptotics of low-degree modes, 
as reflected in the inner phase shift and the small fre- 
quency separations. This is illustrated in Fig. E] for a 
model resembling Procyon. A comparison with the re- 
sults for the much more benign solar model, illustrated 
in Fig. [1] clearly shows the change in the behaviour of 
the inner phase shift Si. This indicates the potential 
for using frequencies of low-degree modes to probe the 
properties of convective cores in stars showing solar-like 
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Fig. 6 The internal phase shift 5i for a model approximat- 
ing Procyon, with mass M = 1.5 Mq and central hydrogen 
abundance X c — 0.05. Continuous curve: I = 0; dotted 
curve: I — 1; dashed curve: I = 2; and dot-dashed curve: 
I — 3. The glitches for / = 2 and 3 are probably associated 
with mixed modes (see Section r2,4|l 



oscillations (see al s o | Popielski fe Dziembowskil 120051 ) . 
Cunha fc Metcalfe! (|2007l ) considered the specific case 
of a small growing convective core, deriving diagnostics, 
based on an asymptotic analysis of the oscillations, of 
the properties of the core. They found that the quantity 



D 



Dm 



(16) 



(see Eqs |4] and E} provides a measure of stellar age; 
more specifically, it is sensitive to the discontinuity in 
the sound speed resulting from the composition discon- 
tinuity at the edge of the core. 

2.3 Diagnostics of acoustic glitches 

The asymptotic expression ([TJ is valid only if the spa- 
tial scale of variation of the equilibrium state is every- 
where much greater than the inverse vertical wavenum- 
ber of the mode. But that condition is not satisfied 
in the Sun and main-sequence solar-like stars: there is 
small-scale variation associated with ionization of abun- 
dant elements and the near discontinuity in low deriva- 
tives of the density at the base of the convection zone. 
These abrupt variations induce small-amplitude oscil- 
latory components (with respect to frequency) in the 
spacing of the cyclic eigenfrequencies v n \ of seismic os- 
cillation and consequently also in Av n i and 5vni- We 
call such abrupt variations an acoustic glitch. 

A convenient and easily evaluated measure of the 
oscillatory component produced by acoustic glitches is 
the second multiplet-frequency difference with respect 



Table 1 Scismically determined prop erties of acoustic 
glitches in the S un (using BiSON data; iBasu et all 120071 ) 



and in Model S (IChristensen-Dalseaard et al. 19961 ) based 



on the analysis by iHoudek fc Goughl (|2007l ) of low-degree 
modes (/ = 0, 3): 571/71 |n and m are the relative depres- 
sion of the first adiabatic exponent and the acoustic depth 
of the second helium ionization zone, and r c is the acoustic 
depth of the sharp transition from radiative to convective 
heat transport at the base of the convection zone 





— 57i/7i n 


Til (s) 


r c (s) 


Sun 


0.043 


819 


2310 


Model S 


0.045 


815 


2270 



to order n amongst modes of like degree I: 

A 2 v n i = v„-n -2v nl + P n+ it (17) 

(lGoughlll990l ). Any localized region of rapid variation of 
cither the sound speed c or the density scale height, or 
a spatial derivative of them, induces an oscillatory com- 
ponent in A2V (the subscripts nl have been dropped) 
with a 'cyclic frequency' approximately equal to twice 
the acoustic depth 



x dr 



(18) 



I-glitch 



(not to be confused with the acoustic radius, Eq.[5]) of 
the glitch, and with an amplitude which depends on the 
amplitude of the glitch and which decays with v once 
the inverse radial wavenumber of the mode becomes 
comparable with or less than the radial extent of the 
glitch. 

Various approximate formulae for the oscillatory 
components that are associated with the hel ium ioniza 



tion h a ve been suggest e d and u sed, by e.g-.lBasu et al 
1994 1 ): liasu fe Antial d2004l). iMonteiro fe Thompson 



1998L l2005h . and iGoughl (j2002l ). not all of which are 
derived directly from explicit acoustic glitches. Gough 
used an analytic function for modelling the dip in the 
first adiabatic exponent 71 = (91np/91np) s , where p, p 
and s are pressure, density and specific entropy. In 
contrast, Monteiro & Thompson assumed a triangular 
form. Basu et al. have adopted a seismic signature for 
helium ionization that is similar to that arising from 
a single discontinuity; the artificial discontinuities in 
the sound speed and its derivatives that this and the 
triangular representations possess cause the amplitude 
of the oscillatory signal to decay with frequency too 
gradually, although that deficiency may not be imme- 
diately noticeable within the limited frequency range 
in which adequate asteroscismic dat a are or will immi- 
nently be available. More recently IHoudek fc Gougb 
(|20071 ) proposed a seismic diagnostic in which the vari- 
ation of the first adiabatic exponent 71 in the helium 
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ionization zone is represented with a pair of Gaussian 
functions. This correctly results in a decay of the am- 
plitude of the seismic signature with oscillation fre- 
quency that is faster than that which the triangular 
and the single-discontinuity approximations imply, and 
also takes some acc ount of the two ion i zation states of 
helium. Moreover, iHoudek fc Goughl (|2007l ) incorpo- 



rated the acoustic cutoff frequency into the variation of 
the eigenfunction phase with acoustic depth r, thereby 
improving the discrepancy between the scismically in- 
ferred depths of the acoustic glitches and that of a cor- 
responding stellar model. 

Fig. [7] shows the outcome of fitting Houdek & 
Gough's (2007) seismic diagnostic to BiSON data (top 
panel), together with the individual contributions from 
the acoustic glitches (lower panel) . The individual con- 
tributions provide a measure of the properties of the 
acoustic glitches in the background stratification of the 
star. One might hope that the variation of the sound 
speed c induced by helium ionization might enable one 
to determine from the low-degree eigenfrequencies a 
measure that is directly related to, perhaps even al- 
most proportional to, the helium abundance, with lit- 
tle contamination from other properties of the struc- 
ture of the star. Table [1] compares some of the proper- 
ties of the acoustic glitches obtained from the BiSON 
data, and fro m adiabatically computed eigenfreq uencies 
of Model S (IChristensen-Dalsgaard et al.lll99a ). using 



the analysis by Houdek fc Goughf "( 2007t h The values 
between the Sun and Model S differ by less than 3%, 
demonstrating the diagnostic potential of this seismic 
analysis using only low-degree modes. 

2.4 Mixed modes 

We have so far considered only modes that are predomi- 
nantly of acoustic nature; these p modes are expected at 
the relatively high frequencies where the stochastic ex- 
citation is most efficient. However, in evolved stars, be- 
yond the end of central hydrogen burning, internal grav- 
ity waves, or g modes, reach frequencies in the range of 
the stochastically excited modes. Their properties are 
characterized by the buoyancy, or Brunt- Vaisala, fre- 
quency N, determined by 



N = g 



1 dlnP dlnp 



7i 



dr 



dr 



^(V ad -V 
V 



V, 



(19) 



where g is the gravitational acceleration. The last 
equality assumes the ideal gas law; V = dlnT/dlnp, 
Vad is its adiabatic value, and V M = dln/^/dlnp, where 
fi is the mean molecular weight. In evolved stars g be- 
comes very high near the compact core; an additional 
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Fig. 7 Top: Symbols show second differences k.iv n x 
(Eq. ll7p for low-degree modes of BiSON data collected over 
a period of 4752 days (jBasu et al.ll2007h . The solid curve 
is Houdek & Gough's (2007) seismic diagnostic, determined 
from fitting it to the data by least squares. The dashed curve 
is a smooth contribution of t±.%v n i, modelled as a third-order 
polynomial in guided by the asymptotic expression JT|, 
and represents near-surface effects. Bottom: Individual 
contributions of the seismic diagnostic. The solid curve dis- 
plays the contribution of the second ionization stage of the 
helium, and the dotted curve the first ionization stage of 
helium. The dot-dashed curve is the contribution from the 
discontinuity in the second density derivative at the base of 
the convection zone 



contribution comes from which is large in the region 
of varying hydrogen abundance left behind after nuclear 
fusion. As a result, N may reach very high values in 
the central parts of the star. 

As an example we consider a model of the sub- 
giant r] Boo, where for the first time individual fre- 
quencies o f solar-like oscillatio ns in another star were 
identified (jKieldsen et al. 1995). Thi s is ba sed on mod- 

thc model 



( 2003h : 



clling of the star by iDi Mauro et al. 
has a mass of 1.7 Mq and a heavy-element abundance 
Z = 0.04. Figure [5] shows the buoyancy frequency at a 
model age of 2.34 Gyr. It is evident that it reaches fre- 
quencies far higher than the acoustic cut-off frequency, 
around 1150 /xHz, in the atmosphere of the model which 
defines the upper limit to the frequencies of trapped 
acoustic modes. Also shown are the acoustic, or Lamb, 
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Fig. 8 Characteristic frequencies for a model of rj Boo. 
The solid curve shows N/2n where N is the buoyancy fre- 
quency (cf. Eg. I19p . and the dashed and dot-dashed curves 
show Si/2tt for I — 1 and 2, respectively. The heavy hori- 
zontal line marks a ty pical observed frequency of r\ Boo (see 
iDi Mauro et al.ll2003h 



frequencies Si and S2, denned by 



S? = 



l(l + l)c 2 



Age (Gyr) 
2.30 2.31 2.32 2.33 2.34 
1 100 E 



2.35 




6200 6180 6160 6140 6120 6100 
T*n (K) 



(20) 



The properties of a mode of oscillation are controlled 
by the magnitude of the frequency relative to Si and 
N. In the region of the star where the frequency ex- 
ceeds N and Si the mode behaves predominantly as an 
acoustic wave. Where the frequency is below N and Si 
the mode has the character of an internal gravity wave. 
Finally, in the so-called evanescent region, where the 
frequency is between TV and Si, the mode amplitude 
varies exponentially with position, cither increasing or 
decreasing. 

In Fig.[S]the horizontal line marks a typical observed 
frequency in 77 Boo. It is evident that this mode may 
have g-mode behaviour in the core of the model and 
p-mode behaviour in the envelope. The dominant be- 
haviour depends on the variation in the evanescent re- 
gion. If the cigenfunction decreases exponentially with 
r in this region the mode has predominantly the char- 
acter of a g mode, whereas if it increases exponentially, 
the mode has the character of a p mode. The discrimi- 
nation between the two cases depends on the strength 
of the exponential variation and hence on the extent 
of the evanescent region. Thus the separation between 
the two cases is stronger for I = 2 than for I = 1 modes. 

With the evolution on the sub-giant branch the star 
expands, and hence the acoustic frequencies decrease 
approximately as p 1 / 2 . On the other hand, the in- 
creasing central condensation increases the buoyancy 



Fig. 9 Evolution of oscillation properties with age (up- 
per abscissa) and effective temperature (lower abscissa) in 
a model of n Boo. Panel (a) shows frequencies of modes with 
/ = (dashed curves) and / = 1 (solid curves). In panel (b) 
the solid curves show the evolution of the dimensionless in- 
ertia E (cf. Eq. I21|l for the two I = 1 modes marked by 
triangles and squares in panel (a), and of a neighbouring 
mode with I = (dashed curve). The heavy vertical line 
identifies the model providing the best fit to the observed 
frequencies. 



frequency and hence the frequencies of the g modes. 
The combined effect for modes of degree I = and 1, 
in models near the identified stage of evolution of r) 
Boo, is illustrated in Fig. For the radial modes, with 
I = 0, and the predominantly acoustic I = 1 modes 
the frequencies decrease with increasing age and hence 
decreasing effective temperature. However, there is ev- 
idently also a branch of frequencies, corresponding to 
a mode of predominantly g-mode character, increasing 
with age. Where this meets a predominantly acous- 
tic mode, the t wo modes undergo an avoided cross- 
ing ( Osakil 1975| l , where the frequencies approach quite 
closely without actually crossing. An example occurs at 
the vertical line, marking the model identified as best 
fitting the observations. At closest approach, the modes 
have a fully mixed nature, with similar amplitudes (in 



Prospects for asteroseismology 



9 



Age (Gyr) 

2.30 2.31 2.32 2.33 2.34 2.35 





1 100 : 




1000 I 


Hz) 


900 ; 








800 k 




700 r 




600 : 




10 6 : 




1(T 



10" 



10" 




6200 6180 6160 6140 6120 6100 
T M (K) 



Fig. 10 As Fig. H but for I = 2. 

terms of energy density) in the central g-mode region 
and the p-mode envelope. 

The variation in the physical nature of the modes can 
be followed by considering the normalized mode inertia 



E 



IyP\SrfdV 
M\8r phot \ 2 



(21) 



where 8r is the displacement vector, £r p hot is its pho- 
tospheric value, and integration is over the volume V 
of the star. This is illustrated in panel (b) of Fig. [5] 
for two modes with I = 1 undergoing an avoided cross- 
ing, and for a neighbouring radial mode. For the latter, 
a purely acoustic mode, E varies slowly with age. The 
I = 1 modes have E close to the radial modes when they 
behave as acoustic modes, while E is somewhat higher 
where the modes have substantial g-mode character. As 
is clear from the figure the two modes exchange char- 
acter at the avoided crossing; at the closest approach, 
near the vertical line, E is the same in the two modes. 

The corresponding results for I = 2 are shown in 
Fig. 1101 Here the number of g-mode branches is higher. 
Again, these undergo avoided crossings with the acous- 
tic modes, but the minimum separation is so small that 
the avoided nature is barely discernible. Also, the mode 



inertias, shown in panel (b) of the figure, on the g- 
modc branches exceed those of the radial modes by 
almost two orders of magnitude (note the logarithmic 
scale). The reason for this difference in behaviour be- 
tween I = 1 and 2 is the greater extent of the evanes- 
cent region in the latter case (cf. Fig. [5]). This decreases 
the coupling between the g-mode and acoustic regions, 
making the g modes much more strongly trapped in the 
innermost parts of the model, and decreasing the inter- 
val over which the avoided c rossing take place (see also 
Christensen-Dalsgaar dl ll980h . 

Oscillations with g-mode character are much more 
sensitive to the properties of stellar cores than are the 
purely acoustic modes. In particular, Eq. (|19p shows 
that the buoyancy frequency depends directly on the 
gradient in the hydrogen abundance. Thus such mixed 
modes have a very substantial asteroseismic diagnostic 
potential. An important issue is whether the modes can 
be expected to be excited to observable amplitudes. It 
is perhaps not unreasonable to assume that the stochas- 
tic excitation excites the modes to a total energy that 
depends predominantly on the properties of the oscil- 
lations in the near-surface layers, where both the exci- 
tation and the dominant damping take place. In these 
superficial layers the characteristic timescale of the con- 
vection is of the same order of ma gnitude as the pulsa- 
tion period of the oscillations (e.g. Goldreich fc Keelev 



1977t IChristensen-Dalsgaard fc Frandsenl Il983h . This 
leads to a strong coupling between the global oscil- 
lations and the turbulent velocity field and conse- 
quently efficient mode excitatio n by the convection (e.g. 
Balmforthlll992bl ; lHoudekl2006[) . Since the properties of 
the oscillations in this region, for low-degree modes, are 
mainly determined by the frequency, we would there- 
fore expect that the mode energy £ = £(u) is also a 
function of frequency. Using that £ ~ EM(V 2 ) (cf. 
Eq. [21]), where (V 2 ) is the mean square amplitude, we 
might therefore expect 



(V 2 ) 1 ' 2 otE- 1 ' 2 



(22) 



at a given frequency. Based on Figs [51 and 1101 we would 
therefore predict that for I = 1 even modes on the g- 
modc branch would have amplitudes that could be ob- 
servable, whereas for I = 2 modes with some g-mode 
character would only be visible in the narrow regions 
around the avoided crossings, for this modclQ 



1 A more careful analysis of the visibility of the modes in terms of 
the peak height in the power spectra, rather than the root-mcan- 
square velocity, shows a somewhat more complex picture, with 
strong dependence on the m ode lifetime in relation to the dura- 
tion of the observations: sec Fletcher ct al. | j2006l) . lDu"pret et al.l 
j2009h . 
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Fig. 11 Echelle diagram (see text), with a starting fre- 
quency of 33.7 /iHz and a large separation Au = 40.3 /iHz. 
Modes of degree I — 0, 1 and 2 are shown with circles, trian- 
gles and squares, respectively. Ope n symbols show result s 
for the best-fitting model of r\ Boo (|Di Mauro et alj|2003h . 
marked by a vertical line in Figs [9] an d 1101 while filled sym - 
bols are for combined ob servations of iKieldsen et all (|2003l ) 
and lCarrier et al.l (120051 ). 



It is of great interest to compare these model prop- 
erties with the observed frequencies. We base this on a 
re-analy sis by H. Kjeldsen (un publi shed) of the observa - 
tion s bvlKieldsen et all |2003T ) and lCarrier et all |2005h 
(see Aerts et alT 20091 for details). The results are il- 
lustrated in an echelle diagram (see lGrec et al.|[l983 ) in 
Fig. 111! Here the frequency spectrum has been divided 
into slices of a length given by the large separation Av, 
starting at a frequency of 33.7 //Hz. These reduced fre- 
quencies are plotted against the starting frequency of 
each slice, corresponding, in the figure, essentially to 
stacking the slices. According to the asymptotic ex- 
pression for acoustic modes, Eq. ([1]), one would expect 
the result to be nearly vertical sequences of points, sep- 
arated by the small separation, although perhaps with 
some variation arising from the acoustic glitch associ- 
ated with the second helium ionization zone (cf. Sec- 
tion [23]). For the observations, shown by filled sym- 
bols, this is indeed satisfied by the frequencies for I = 
and 2. However, the points for I = 1 clearly show a less 
regular behaviour. 

To interpret these results the open symbols show cor- 
responding results for a model matched to the observed 
frequencies, as well as the effective temperature and 
luminosity of the star; this model corresponds to the 
vertical lines in Figs H] and [lJJ] The frequencies were 
corrected for near-surface effect s according to the pro- 
cedure of iKjeldsen et al.l (|2008h (see Section |2~5]) . To 
provide some indication of the visibility of the modes, 
the size of the symbols has been scaled to reflect the am- 
plitude estimate in Eq. (|22|) . relative to a radial mode 



of the same frequency. For I = 2 this essentially makes 
the modes on the g-mode branches invisible, owing to 
their large E. However, the I = 1 modes are retained, 
with a behaviour showing strong departures from the 
acoustic asymptotic relation, for modes with a strong 
mixed nature. As an example, the dotted line connects 
the two modes shown to undergo an avoided crossing 
in Fig. as expected for modes at the closest point of 
an avoided crossing these are distributed equally on ei- 
ther side of the expected location of the acoustic mode. 
Although there is no detailed agreement between the 
observed and computed frequencies of the I = 1 modes, 
the strong resemblance between their behaviour sug- 
gests that the observations show evidence for mixed 
modes in this star. 

2.5 Superficial problems 

A major issue in the analysis of solar- like oscillation fre- 
quencies is the influence of the near-surface layers. Stel- 
lar modelling normally treats these layers in a highly 
simplified manner, with a simple model of the stellar 
atmosphere and a simplified treatment of the upper, su- 
peradiabatic, parts of the c onvection zone, such as the 
mixing-length formulation ( Bohm-Vitensd 1958 1. Also, 
the effects of the turbulent Reynolds stresses on the 
structure are almost always ignored. In computations 
of oscillation frequencies the adiabatic approximation is 
commonly adopted, and perturbations to the turbulent 
Reynolds stresses are neglected. When nonadiabatic 
calculations are carried out the treatment of the per- 
turbation to the connect ive flux is a major uncertainty 



in the calculation (e. 
198dlBalmforthlll992 



Baker fc Goughl Il979l : iGoueh 



These problems arc confined to the stellar atmo- 
sphere and the outermost part of the the stellar in- 
terior. In these regions the properties of modes of 
low or intermediate degree are largely independent of 
degree, at least in spherically symmetric models, and 
hence the influence of the near-surface problems on 
the oscillation fre quency is essentially just a function 
of fre quency (e.g. JChristensen-Dalsgaard fc Thompson 
19971 ): for modes of intermediate degree there is some 
dependence on degree in the form of a trivial scaling 
with mode inertia. In the solar case, where modes over 
a broad range of degree are observed, these properties 
allow the near-surface errors in the model to be isolated 
and suppressed in the analysis of the observed frequen- 
cies. In fact, for models such as Model S they dominate 
the d ifferences between the solar and mo del frequencies 
(e.g.. Christensen-Dalsgaard et al.lll996 ). This provides 
an estimate of the effects in the solar case, showing 
that they are very small at low frequency and increase 
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rapidly with fr equency, as indeed expected from simpl e 
analyses (e.g., IChristensen-Dalsgaard fc Cxoughlllism . 
As a consequence, they also have a substantial effect 
on the determination of the large frequency separation 
(see Eq.i]). 

Such a separation is not possible for stellar observa- 
tions, restricted to low degree; hence the near-surface 
effects must be kept in mind as a source of systematic 
error in the analysis. As discussed in Section [221 the 
effects are strongly suppressed in separation ratios (cf. 
Eq. [TU, which are essentially determined just by the 
inner phase shift. This is extremely valuable for the di- 
agnostics of stellar cores. To obtain further information 
from the observed frequencies, additional assumptions, 
so far with limited justification, must be made. The 
behaviour of the oscillations in the stellar atmosphere 
is largely controlled by the value of the frequency, mea- 
sured in terms of the atmospheric acoustical cut-off fre- 
quency w ac which, in the approximation of an isother- 
mal atmosphere, is given by 



c 
2H 



(23) 



where H is the pressure scale height (note that in an 
isothermal atmosphere with constant mean molecular 
weight the pressure scale height is equal to the den- 
sity scale height). Thus, for stars similar to the Sun 
one can perhaps use the ansatz that the near-surface 
effects can be approximated by afjQ (w/tj ac ) where 
/surf * s obtained from the solar observations and a is 
a constant to be determined as part of the fit of the 
model to the observe d frequencies. In a similar spirit 
Kjeldsen et al. ( 2008| ) suggested to represent the near- 
surface frequency shift as a(w/wo) fc , where lu is a suit- 
able reference frequency, the exponent b is obtained 
from the solar data, and a is obtained from the fre- 
quency fit, together with a scaling factor that accounts 
for the difference in mean density between the star and 
the model. 

It is obvious that these simple procedures have rela- 
tively weak justification, particularly for stars that are 
somewhat different from the Sun. Analysis of the de- 
tailed data for a broad range of stars that are being ob- 
tained by the Kepler mission will certainly contribute 
to a better characterization of the near-surface effects. 
Also, an overview of how these effects influence the fre- 
quencies of stellar models, depending on the stellar pa- 
rameters, would be very interesting. 

A key goal is clearly to improve the treatment of 
the physics of the near-surface regions. In the solar 
case, models of the near-surface region based on hy- 
drodynamical simulations result in frequencies in bet- 
ter, but still not complete, agreement with th e obser- 
vations dRose^haTeTaD [1993; [lEsT^ EqqI). Such 



calculations were also applied with apparent success 
to the analysis of MOST observations of r] Boo by 
Straka et alj ((2006) . The effects of convection dynam- 
ics on the oscillati o n frequencies were also disc ussed by 
Balmforthl |l992bl ). lRosenthal et al1(|l995l ) and lrloudek 
( 19961 ). who used the nonlo cal, time-dependent for- 
mulation by iGoughl (|l977al |b1) in their linear stabil- 
ity analyses. These authors reported that the in- 
clusion of the Reynolds stress in an equilibrium so- 
lar model reduces substantially the frequency differ- 
ence between observed and adiabatically computed ra- 
dial eigenfrequencies. This difference, however, be- 
comes larger again when nonadiabaticity and convec- 
tion dynamic effects are included in the stability com- 
putations, thereby compensating nearly the frequency 
shifts arising from the Reynolds stress in the cquilib- 
ri um model. A similar conclusio n was reported by 



Christensen-Dalsgaard et al.l (|1995l ) for r\ Boo. The ef- 
fect of such modelling on other stars should obviously 
be investigated. 



3 The observational situation 

The tiny amplitudes of solar-like oscillations have 
greatly complicated their observation in distant stars, 
despite extensive effor ts. It was only the observa- 



tions of iKieldsen et alJ (|1995l ) that led to the first de- 
tection of individual modes of oscillations in a solar- 
like sta£ J __the_star 77 Bo o; this was later confirme d 



Carrier et all j2005f) 



by IKieldsen et all (|2003l ) and 

However, in the last decade the development in ob- 
servational techniques to measure Dopplcr velocity, 
largely motivated by the search for extra-solar plan- 
etary systems, has led to fairly detailed observations 
of solar-like oscill ations in a number of sta rs (for a 



brief review, see Bedding fc Kjeldsen 2008h . reach 



ing a noise lev el in the ampl i tude spectrum as low 



as 



1.4 cms 1 (jKieldsen et al.l 120051 ). From space, 



photometric observations allow simultaneous observa- 
tions of large numb ers of stars. The C0R0T mission 
jBaglin et all l2009h las produced remar kable astero- 



seismic results (e.g.. iMichel et al.l l2008al lrJ). including 



striki ng results on solar-like pulsations in red giants 



(e.g., iDe Ridder et~all 120091). Also, the firs t results 



from the Kepler mission (jBorucki etaUl2009l ) confirm 
the g reat potential of the mission for as teroseismology 
fe.g.. lChristensen-Dalsgaard et al.ll2007l h 

The observations generally provide a reliable mea- 
sure of the large frequency separation (see Eq. fj]) . To 
move beyond this an identification of at least the degree 
of the modes is required. In the ca se, for example, of 
the binary system a Cen A and B (jBouchv fc Carrier 
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Bazot et all2007tlCarrier fc BourbaiJ200d|Kieldsen et al 



20011 l200a Butler et al. 200 4; Bedding et all 12004 



200E 



^) or for (3 Hydri ( Bedding et ail 2007 ) an unam 



biguous identification has been possible, based on a 
match to the asymptotic behaviour (cf. Eq. [I]). This 
has led to detailed asteroseismic analyses for these 
stars, yielding precise determinations of their over- 
all parameters and some indi cations of possible prob 



Eggenberger et al.l 12004 : 



lems in the modelling (e .g 
Miglio fc Montalbanll200j ). 

In other cases, however, the mode identification has 
proved to be difficult. An interesting example is the star 
Procyon A (a CMi), of spectral type F5IV-V, which 
has been the target of several asteroseismic observing 
campaigns. Evidence for a power excess as sociated with 



solar-l ike oscillations was found already bv lBrown et al 
1 19911 ): however, several subsequent observations have 



failed to make an unambiguous identification of the 
modes of oscillation, and doubt was even raised on the 
reality of the oscillations on the basis of observ ations 
from the MOST satellite jMatthews et al.ll2004l ) (see, 



however, Bedding et al. 2005h . A very extensive ob 



serving campaign, inv o lving 1 1 telescopes, was carried 
out by lArentoft et alJ (|2008l ) . This confirmed beyond 
any doubt that Procyon exhibited solar- like oscillations; 
the mode identification has proven to be difficult but is 
now close to being completed (Bedding et al., in prepa- 
ration). Similar problems have been found in the anal- 
ysis of data from the C0R0T satell ite, such as the star 



HD49933, of spectral type F2V (lAppourchaux et al 
2008HBenomar et ahlbooa iGaulme et al.ll2009h ~ 



These problems in the data analysis appear to a large 
extent to be associated with the mode lifetime in these 
stars, which is substantially shorter than in, for exam- 
ple, the Sun or a Cen A and B, hence leading to a more 
complex spectrum of oscillations. Thus it is of substan- 
tial interest that a close relation has been found be- 
tween the effective temperature T e ff of low-mass main- 
sequences stars and the average mode lifetime (r) of 
the five most prominent modes: 



(24) 



(jChaplin et alJ 120091 ). This relation can obviously be 



of importance in the simulation of asteroseismic oscilla- 
tions and the selection of targets. However, it should in 
any case not be a surprise that further development of 
the analysis techniques is required for stars whose prop- 
erties differ substantially from those of the Sun. These 
developments will of course be inspired by the exten- 
sive data obtained from C0R0T and expected from Ke- 
pler. However, analysis of realistically simulated data 
for a variety of stellar parameters is also very valu- 
able, ideally carried out in blind, so-called 'Hare and 
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Fig. 12 Relative differences in squared sound speed be- 
tween the Sun and three solar models, in the sense (Sun) - 
(model), inferred from inver sion of the so-ca lled 'Best set' 
of observed frequencies (see iBasu et al.l Il997l ) . The barely 
visible vertical bars mark 1 a errors in the inferences, while 
the horizontal bars indicate t he resolution of the inversion. 
Open symbols: Model S of IChristensen-Dalsgaard et al.l 
dl996|); filled circl es: corresponding model, but assuming 



the lAspIundl (|2005l ) compositi on; dashed curve: cor respond- 
ing model, but assuming the lAsplund et all (|2009l ) compo- 
sition; dot-dashed curve: model corresponding to Model S, 
but switching off electron screening (in the core this can 
hardly be distinguished from the dashed curve) 



Hounds', investigatio ns. Tools have been dev eloped for 
such simulations (e.g., De Ridder et al.|[2006l ). and they 
are being used in extensive analyses in the asteroFLAG 
project, with specific ap plication to the data expected 
from the Kepler mission (jChaplin et al.l 120081 ) . 



4 Effects of model physics 

An ultimate goal of helio- and asteroseismology is ob- 
viously to relate the inferences to the physics of stellar 
interiors and hence to get a better understanding of 
the processes taking place in stars. Here we consider a 
few relevant examples, including some that were high- 
lighted at the workshop, and consider effects on stellar 
structure that arc, or may be, amenable to helio- and 
asteroseismic investigations. 

4.1 The solar sound speed 

The observed solar oscillation frequencies cover a broad 
range in spherical-harmonic degree. This has allowed 
inverse analyses to determine solar internal properties, 
in particular the sound speed, with high accuracy in 
most of the solar interior. Solar models computed a 
decade ago generally provided relatively good agree- 
ment with these inferences, with maximum relative de- 
viations in the square of the sound speed well below 
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Fig. 13 Changes to the natural logarithm In « of the opac- 
ity, required to bring models with revised composition in 
accord ance with Model S of IChristensen-Dalseaard et al.l 



l| 19961 ). Solid curve: correction for the 
Dosition; dashed curve: correction for the 
(120091 ) composition. 



AsplundlT2005l) com 



Asplund et al 



(See IChristensen-Dalseaard et al 



1 per cent. A typical example is illustrated by the 
open symbols in Fig. [T2l based on the so -called Model 
S of IChristensen-Dalsgaard et al. ( 1996 ) ; in common 
with many solar models of the epoch it assumed the 
Grevesse fc Noelsl (Il993h solar compositi on. However, 
starting with lAllende Prieto et al.l (|200l[ ) reanalyses of 
solar spectral lines have led to substantial revisions in 
the determination of the surface abundances of, in par- 
ticular, carbon, nitrogen and oxygen, leading to a re- 
duction in the ratio Z s /X s between the abundances by 
mass Z s of heavy elements and X s of hydrogen from 
0.0245 to 0.0165. This revision resulted from the use of 
hydrodynamical models of the solar atmosphere and the 
inc lusion of NLTE effects in the analy sis (for reviews, 



Asplundlliool lAsplund et al.ll2005l ). The main ef- 



scc 

feet on solar modelling is a corresponding reduction in 
the opacity, leading to a substantial change in the sound 
speed in the radiative interior. The effect on the sound 
speed is illustrated by the filled symbols in Fig. mil it 
is evident that this has led to a serious deterioration in 
the agreement between the solar models and the helio- 
scismic inference. 

This conflict between helioseismology and the solar 
models resulting fr om the revised comp osition was dis- 
cussed in detail bv lBasu fe Antia ( 20081 ); possible solu- 
tions were reviewed by Guzikl ( 2006l k with the conclu- 
sion that none of the proposals were fully satisfactory. 



2 Here, and in the following results on solar modelling, the mod- 
els correspond to Model S except in the specific property under 
investigation; the models have all been calibrated to solar ra- 
dius, luminosity and the assumed value of Z B /X B at solar age, by 
adjusting the mixing length and the initial composition. 



Since the dominant effect on the models of the change 
in the composition arises from the resulting change in 
opacity, a simple solution would be to ma ke correspond 



ing c ha nges to the intrinsic opacity ( e.g., Bahcall et al 



20051 ). Christensen-Dalsgaard et al. ( 2009 ) estimated 



the opacity change required to restore the model com- 
puted with the old composition; the result is illustrated 
in Fig. 1131 It is perhaps questionable if such a major 
change, of up to 30 per cent, is physically justifiable, al- 
though there are undoubtedly significant uncertainties 
in the opacit y calculations. 

Recently, I Asplund et al.l (|2009l ) reconsidered the 



abundance determinations, using improved solar atmo- 
sphere models and a car eful selection of atomic data and 
spectral lines (see also Grevesse et al. 2010l this vol- 
ume) . This resulted in a modest increase in the inferred 
abundances, such that Z s /X s = 0.0181. The effect on 
the comparison between the resulting solar model and 
helioseismic inferences is illustrated by the dashed line 
in Fig. [121 Obviously the model is somewhat closer to 
the solar sound speed although the differences are still 
much larger than for the original model. The opacity 
change required to restore the original model is illus- 
trated in Fig. 1131 the maximal change, around 23 per 
cent, may still be larger than can reasonably be ex- 
plained by errors in the opacity determinations. 

Other uncertainties in stellar modelling may have a 
noticeable effect on the comparison with helioseismic 
inferences. An interesting example is the effect of elec- 
tron screening on the nuclear reactions. This is nor- 
mally treated in the Debye-Hiickel ap proximat i on, ac - 
cording to a formulation developed bv lSalpeterl |l95J); 
here the formation of clouds of negative charge around 
the reacting nuclei leads to a reduction in the Coulomb 
repulsion and hence an increase in the reaction rates. 
However, it has been pointed out that the underlying 
mean-field approximation may be invalid un der stel- 
lar conditions (see, for example, Shaviv 20041 for a re- 
view). As discuss ed by iMussack fe Dappe X ifeoiol . this 
volume) (see also iMao et al.ll2009l) this has motivated 
calculations of the properties of stellar plasmas using 
the techniques of molecular dynamics. The results in- 
dicate that the effects of electron screening are weaker 
than in the Salpeter formulation and that the interac- 
tions between the charged particles may in fact lead to 
a reduction in the reaction rates. To investigate the 
consequences for solar modelling we have computed a 
model where electron screening was simply switched off. 
The resulting sound-speed difference is illustrated by 
the dot-dashed line in Fig. [TSJ Although less dramatic 
than the effect of the change in composition, the change 
in the nuclear reactions, and the resulting change in 
the composition required to calibrate the model, leads 
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to a definite deterioration in the agreement with the 
helioscismically inferred sound speed; such an effect of 
su ppressing ele c tron s creening was already pointed out 
bv I Weiss et ail |200lh . 



4.2 The onset of convective cores 

A major uncertainty in the modelling of intermediate- 
and high-mass stars is the treatment of convective 
cores. Convective cores are typically found in stars of 
masses just a little higher than the mass of the Sun, thus 
including also stars where solar-like oscillations are ex- 
pected. As discussed in Section |2"T2"1 the detailed proper- 
ties of the core are reflected in the small frequency sep- 
arations and hence there are excellent prospects to use 
asteroseismology to constrain the modelling of convec- 
tiv e cores. A detailed analys i s of th is was carried out 
by Popielski fc DziembowskJ ( 2005h . emphasizing the 
effect of the composition discontinuity caused by the 
growing convective core in stars of mass up to around 
1.9 M Q . 

Here we consider the effects of modifications to the 
microphysics and composition of the stellar interiors. 
The convective instability is caused by the increasing 
temperature sensitivity of nuclear reactions as the CNO 
cycle, with increasing stellar mass, becomes dominant 
in hydrogen burning. An additional contribution is the 
gradual conversion of 16 O to 14 N through a branch of 
the CNO cycle; this increases the abundance of 14 N 
which controls the overall rate of the energy generation 
from the CNO cycle and hence its importance, leading 
to the growth of the mass of the convective core. 

It is obvious that the CNO cycle is sensitive to the 
abundances of carbon, nitrogen and oxygen. Abun- 
dances of solar-like stars are typically determined rela- 
tive to the solar abundance, and hence the recent revi- 
sions of the solar abundances have a similar effect on the 
assumed abundances of stars and h ence potentially on 



the onset of convective cores. In fact, VandenBerg et al 



(2007J) showed that the lAsplund et all (|2005h abun 



dances led to isochrones for the open cluster M67 which 
did not show a 'hook', resulting from the presence of a 
convective core. This was in contrast to models with 
the old composition, which showed a clear hook, and in 
apparent contradiction to observed colour-magnitude 
diagrams for the cluster. To investigate this further, 
we have considered the properties of convective cores 
as a function of stellar mass, characterizing them by 
the maximum mass reached by the core during central 
hydrogen burning. The models neglected diffusion and 
settling and used the NACRE (jAngulo et al.lll999h nu- 
clear param eters, but otherwise largely co r respon ded to 
Model S of IChristensen-Dalsgaard et all (|l996h . The 




Fig. 14 Maximum fractional mass of the convective core 
during main-sequence evolution, as a function stellar mass in 
units of the solar mass. The models used the lAngulo et al.l 
1 19991 ) nuclear parameters and included electron screening; 
diffusion and settlin g was not included. Sol id line: refer- 
ence model, with the lGrevesse fc Noels! (|l993l ) composition, 
and initial hydrogen and heavy-element a bundance s given 
by (X ,Z ) = (0.708, 0.0196); dashed line: lAsplund! (|2005h 
comp osition, and (Xq,Zq) = (0.730,0.0137), dot-dashed 
line: lAsplund et alj ( 20091 ) composition, and (Xo,Zo) = 
(0.722,0.0149) 



results are illustrated in Fig. 1141 In agreement with 
VandenBerg et alJ (|2007l ) the mass at onset of the con- 
vec tive core increa s es by around 0.04 M Q as a result of 
the lAsplund et alJ (l2005|) revision of the com position. 
The modest revision by lAsplund et al.l (|2009l ) has only 
a slight effect, in comparison. 

The presence of convective cores is obviously also 
sensitive to nuclear properties, particularly for the reac- 



tion 14 N( 1 H, 7) 15 which controls the rate of the dom- 
inant part of the CNO c ycles. From new la boratory 



experiments and analyses. I Angulo et alJ (|2005l ) found a 
large reduction in this rate. As shown in Fig. [15] this 
has a substantial effect on the minimum mass of stars 
with convective cores, increasing it by 0.06 M Q . On the 
other hand, since the CNO cycle makes a modest con- 
tribution to solar energy generation the effect on solar 
models is largely insignificant. 

We have also considered the effect on the convective 
cores of switching off electron screening; the effect of 
electron screening increases with the charges of the in- 
teracting nuclei and hence turning it off decreases the 
relative importance of the CNO cycle, again delaying 
the onset of convective cores. This is also illustrated in 
Fig- El the effect is relatively modest. 

Further investigations will be required to test the 
effects on the oscillation frequencies of these changes 
in the model properties, compared with other uncer- 
tainties in the treatment of convective cores, and the 
prospects for applying asteroseismic diagnostics, such 
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Fig. 15 As Fig. ITU all models having the 
iGrevesse fc Noels! (|l99ot ) composition. Solid line: ref- 
erence model; das hed line: reduced 14 N reaction rate 
(|Angulo et alj|2005h : dot-dashed line: no electron screening 



as those proposed by Cunha fc Metcalfe! ( 2007 ). to in- 
vestigate them. 



an example is the Stellar Observations Network Group 
(SONG) project (jOrundahl et al.ll2009h . aiming to es- 
tablish a network of 8 observatories, with a suitable 
world-wide distribution, equipped with 1 m telescopes 
and high-resolution spectrographs to carry out astero- 
seismic Doppler-velocity observations. The proposal to 
set up the SIAMOIS facility for Do ppler-velocity obser - 
vations on Dome C in Antarctica (jMosser et alJ 120081 ) 
is also very promising. 

In parallel with these great observational strides 
there is a strong need for further development of the 
techniques for stellar modelling and the interpretation 
of the asteroscismic data. It is already clear that, not 
surprisingly, the oscillation properties of other stars 
show striking departures from the familiar properties 
of solar oscillations and that new techniques will be 
needed to optimize the analysis of the data and the use 
of the results for asteroscismic inferences. From the 
results of these analyses we can surely expect further 
surprises and identification of deficiencies in our stellar 
modelling, leading to new insights into the physics of 
stellar interiors. 



5 Prospects for the future 

Asteroseismology is being revolutionized in the present 
period. As documented in a recently published special 
Volume 506 of Astronomy and Astrophysics, the initial 
results from the CoRoT mission are very promising, 
and at the time of writing the very early data from the 
Kepler mission are being analysed, through an intense 
and coordinated effort in the Kepler Asteroscismic Sci- 
ence Consortium, with a view towards publication in 
early 2010. Given the continuing operations of CoRoT 
and the very extensive data to be obtained with Kepler 
it is clear that coming years will see a huge effort in 
the analysis of asteroscismic data, with a correspond- 
ing increase in our information on and, one may hope, 
understanding of stellar properties. On a somewhat 
longer timesc ale there is hope that the PLATO mission 
1 Catalal 120081 ) will provide detailed asteroseismic data 
on tens of thousands of stars. 

Despite the impressive success of space-based aster- 
oseismic investigations, there is still a need for ground- 
based obser vations . As found early in observations 
of the Sun (Harvey 1988h the intrinsic stellar noise is 
much higher, relative to the oscillations, in photomet- 
ric observations than in observations of Dopplcr veloc- 
ity. Velocity observations can be carried out from the 
ground, the main limiting factor being the access to 
the required instrumentation, at several sites to ensure 
continuous data and over a sufficient length of time. 
This motivates the establishment of dedicated facilities; 
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